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We study the collisional processes that can lead to thermalization in one-dimensional systems. 
For two body collisions excitations of transverse modes are the prerequisite for energy exchange 
and thermalzation. At very low temperatures excitations of transverse modes are exponentially 
suppressed, thermalization by two body collisions stops and the system should become integrable. 
In quantum mechanics virtual excitations of higher radial modes are possible. These virtually 
excited radial modes give rise to effective three-body velocity-changing collisions which lead to 
thermalization. We show that these three-body elastic interactions are suppressed by pairwise 
quantum correlations when approaching the strongly correlated regime. If the relative momentum 
k is small compared to the two-body coupling constant c the three-particle scattering state is 
suppressed by a factor of (fc/c) 12 , which is proportional to 7 12 , that is to the square of the 
three-body correlation function at zero distance in the limit of the Lieb-Liniger parameter 7 > 1. 
This demonstrates that in one dimensional quantum systems it is not the freeze-out of two 
body collisions but the strong quantum correlations which ensures absence of thermalization on 
experimentally relevant time scales. 
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I. INTRODUCTION 

One-dimensional (ID) systems [H,[2| are a model to study the fundamental processes of dynamics and (de)coherence 
in interacting many-body quantum systems. Ultracold atoms in strongly elongated traps with u> r ^> lu z (w r , uj z 
being the frequencies of the radial and longitudinal confinement, respectively) offer the possibility to implement ID 
quantum physics if both the temperature T and chemical potential fi are small compared to the energy scale given 
by the transverse confinement: 

/i < fiuj r , ksT < hui r . (1) 

ID systems of ultra-cold atoms were implemented in both optical lattices Q and atom chips In the limit of zero 
temperature they are a realization of the Lieb-Liniger model [5] of spinless bosons with contact (point-like) interaction, 
a prime example of an integrable system. 

An important parameter characterizing an ID system of bosons with point-like interactions described by the (three 
dimensional) s-wave scattering length a s is the Lieb-Liniger parameter [j| 

7=-^, (2) 

where m is the mass ot the bosonic atom, n\D the linear density of the atoms in the ID trap characterized by the 
transverse confinement frequency uj t , and l r is the fundamental length scale of the localization of an atom in the 
transversal direction given by 

l r = \Jt\j (mu) r ). (3) 

The limit 7^1 corresponds to a weakly-interacting regime, whereas 7 > 1 signifies strongly-interacting, strongly 
correlated (Tonks-Girardeau) regime. 

In an integrable system Jy, 0] the number of their integrals of motion equals exactly the number of their degrees 
of freedom. Thus such a system always "remembers" its initial state in the course of its dynamical evolution, and 
thermalization does not occur. In an integrable system the finite spread of initial energy may lead only to relaxation 
towards the generalized Gibbs (or fully constrained thermodynamic) ensemble Q. Strictly speaking, there is no 
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thermalization in any closed system, but for non-integrable systems the eigenstate thermalization hypothesis Q 
holds, enabling dephasing to mimic the relaxation to the thermal equilibrium. 

Recently strong inhibition of thermalization, was observed on an optical lattice experiment with bosons deep in 
the ID regime [l0|. On the other hand interference experiments on atom chips with pairs of weakly interacting Bose 
gases fulfilling the conditions of Eq. ([I} are in a good agreement with the thermal-equilibrium description of the ID 
atomic ensembles 

In the present paper we investigate the collisional properties of Bose gases in a ID geometry and how they contribute 
to thermalization. We follow thereby a procedure outlined in our two previous works [lil Il5j |. and give a more in- 
depth description of the underlying theoretical considerations. We start with the calcualtion of the freeze-out of 
thermalization providing two body collisions. We then proceed to show that virtual excitations of excited states 
lead to effective three-body collisions, which lead to a term in the Hamiltonian that breaks integrability and enables 
thermalization. We then proced to estimte the effects of quantum correlations in ID and show how they suppress the 
three-body term for strongly correlated ID systems, thereby extending the time scale, on which a quasi-lD system 
can be approximately described as integrable. 

In our theoretical considerations we consider identical bosons in a tight ID wave guide with radial confinement 
given by a 2D harmonic oscillator with frequency w r (we set uj z =0). The contact interaction is described by the 

). The Hamiltonian of the ID system is: 
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Thereby the field operators ip{ r ) are assumed to vanish for 
L. For the solutions to Eq. [S]we make the Ansatz: 
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where L is the quantization length and the atomic annihilation and creation operators a^ n ^ and a| n g, k obey 

the standard bosonic commutation rules. (f n ,i{x,y) is the normalized eigenfunction of both the radial confinement 
Hamiltonian, 

H (r) (pnA x ^y) = {n + l)hu r ip n x(x,y) 

and the z-projection of the orbital momentum, 

~i[x(d/dy) - y(d/dx)](p n ,e(x,y) = £(p n . e (x,y). 

Because we consider identical Bosons the main quantum number n = 0, 1, 2, and the z-projection quantum 
number £ of the orbital-momentum is restricted by 



\£\ = mod (n, 2), mod (n, 2) + 2, 
and thus has the same parity as the main quantum number. 



n — 2, n 



II. TWO-BODY COLLISIONS 



We first look at collisions of two identical bosonic atoms that are initially in the transverse ground state of the 
radial confinement. If the collision is restricted to ID, that is both atoms remain after the collision in the transverse 
ground state, then there can be no energy exchange and cinsequently no thermalization. For such two-body collisions 
to contribute to energy exchange and thermalization, they have to lead to a change in transverse excitation. By 
symmetry Ani + An2 must be even. For atoms in the transverse ground state {n% = ri2 = 0) Ani = A712 band 
following the above considerations their orbital-momentum quantum numbers after collision are restricted to — £ and 





FIG. 1: Feynman diagrams for the effective three-body processes in the second (a) and third (b) orders of perturbation theory. 
Solid and dashed lines correspond to atoms in the ground and excited states of the radial trapping Hamiltonian, respectively. 



+£. The rate of populating the radially excited modes by pairwise atomic collisions T2b, can then be estimated for a 
non-degenerate Bose gas, using Fermi's golden rule. For ksT < fiw r this rate is 



r 



2y/2h 



26 



n 1D a; _• 
e 



mil 



(7) 



The dimensionless quantity 



C = ni D a 2 Jl r 



combines two dimensionless parameters characterizing a ID system. nuja s cx 



7i uj r 



is a measure how much the 



interaction energy (the chemical potential fi) is smaller than the energy scale given by the transverse confinement. 
a s /l r characterizes the relation between the transverse confinement and the strength of the contact interaction. Its 
importance can be seen when looking at how the effective ID coupling constant g±o of pairwise interacting bosonic 
atoms in a waveguide changes with confinement due to virtual excitation of the radial modes. Following Olshanii [l6| 
2huj r a s /[l — c ^ l 3 ], C sa 1.46 and increases as a s /l r grows. This points to the ratio a s /l r < 1 as the measure 
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how much the ID approximation is violated. In a general case we have to change in all the following expressions a s 
to^/Il-^]. 

Eq. (O has a also transparent physical interpretation: The rate is related to the 3D atomic density (~ nw/l 2 ), 
times the s-wave scattering cross-section (~ oq), times the exponential Boltzmann factor for the fraction of atoms 
fast enough to scatter into higher radial modes, times the corresponding velocity of the collision (~ Ti/{ml r )). 

Looking at the scaling of Eq. one immedeately sees that the collision rate leading to thermalization (T2b) 
rapidly diminishes when the temperature approaches T ~ fuxi r and is suppressed by more then a factor of 50 (e~ 4 ) 



Estimating the numbers for recent experients [12| in Rb: a s = 5.3 ran, nm = 50 /xm , cj r /(27r) = 3 



for T — ^huj r 

kHz, T = 30 nK (( ps 0.007) one obtains a collision rate of T2b ~ 0.02 s . Ths is at least one order of magnitude 
too small for two body collisions to be responsible for the thermalization required in the evaporative cooling process 
leading to these low temperatures. 



III. THREE-BODY COLLISIONS 



If the kinetic energy of the collision is less than 2fku r , then the radial modes can be excited only virtually. Such 
processes contribute to the system dynamics in the second and higher orders of perturbation theory. 

The simplest case is when after the collision the radial motion state is |{n'i,^i}, {n-2^'2}) = 1^0, 0}, {2p, 0}). Then 
only one more collision is enough to de-excite the radial mode and bring the system back on the energy shell [see Fig. 
ffl^a)]. Such a process yields an effective three-body collision already in the second order of perturbation theory. 

In contrast processes involving a virtual excitation to £}, {n' 2 ,+£}), £ j£ 0, shown in Fig. QJb), contribute 

only in the third order, and thus will be neglected. 
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A. Perturbative approach 

We will now calculate the matrix elements for the process shown in figure [lja) which leads to effective three-body 
collisions. In our perturbation calculation the small parameter is nuja s , that is the mean field interaction /_i is much 
smaller then the energy scale Tikj r connected to the transverse confinement. In addition, to avoid complications related 
to the confinement-induced resonance in ID scattering [l6j we assume a s <C l r . We can then rewrite the Hamiltonian 
© as 

fh 2 k 2 \ t 

HaD = <L [^r+ nhu} r)a {n e}k a {n , e}k + 



n.l.k 
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where all the terms irrelevant to the process under discussion (Fig. la) are gathered in TZ, and the radial matrix 
element /o^o.o' is S iven b y ; 

/o,o;o, o° = J dx J dy J dx ' J d y'vn=2v4=a( x 'y)vtA x 'iy') 5 ( x ~ x ') 5 (y-y') x 

vqA x ' >y')vofl(. x iy)i p = o, i, 2, 3, ... (9) 

It connects to two atoms in the ground state of the incoming channel, to one atom remaining in the same state, and 
the other being excited to a state with zero orbital-momentum quantum number and even main quantum number 
ri = 2p, p = 0, 1, 2, ... (remember for Boson n and I are required to have the same parity). 

To evaluate Eq. ©, we recall that the normalized radial wave functions Lp n ^{x, y) of interest are real and can be 
expressed through Laguerre polynomials L p 



;„-i„,-«l>:!» (O^exp ( -^±f) Lp ( V*LEZ ) ( J()) 



Since L p (0) = 1, we obtain [I?} 

5{x- x')5(y -y')(po t o{x' ) y')(po ! o(x ) y) = (11) 
independently of p. Then we easily obtain the necessary matrix element as /q o ; o ° = 

;° ; °;°/27r?2 j where the 

coefficient C^p'o'o'o is defined by the expansion 

' x — x 1 y — y'\ fx + x 1 y + y n 

)=■ 

n i 



- x V-y \ x-rx yjry \ ^ />2p,o ; o.o / x r > >\ 

) {~J2~'~JT~ ) = l^l^^p-nA^-e^P-^'yWn^eix ,V )■ (12) 



Comparing the coefficients in front of (x 2 + y 2 ) p in the left- and right-hand sides of Eq. (fl~2|) . we obtain C^p'o'o'o = 2 p 
and 

f 2p,0;0,0 _ 1 n o\ 
/0,0;0,0 - 2P+T-nl 2 ' 

In our consideration we are only interested in the case where the collision energy of the two atoms is always much 

smaller than huj r . Then, using the matrix element (|13j) and adiabatically eliminating the radially excited mode 
operators, we obtain from the original Hamiltonian Eq. (|4]) an effective ID Hamiltonian: 

» v-^ 1i 2 k 2 „+ A huj r a s \ „+ A + - - £,h(jj r a 2 -t -t »t - - - i-x a-k 

UlD = <L ~2m~ a k ak + - — 2w a k+ q a k'- q a k' a k 1^ a k [ a k > 2 a k ' 3 a ^ a ^ a ^> (14) 

k k,k' ,g {kj'} 
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where we write for simplicity = <X{o,o} fc an d the numerical constant £ is given by 

oo 

£ = 4^1/(4» = 41n(4/3)«1.15. (15) 

P =i 

Note that the relative contribution of the virtual states with the excitation energy higher than 2huj r given by (£ — l)/£ 
is remarkably small. The summation in the last term of Eq. (|14j) is taken over all the kinetic momenta obeying the 
conservation law 

K + k' 2 + k' 3 = ki + k 2 + k 3 . 
Introducing the field operator ^p(z) = L^ 1 / 2 ^ k &k exp(ifcz), we can rewrite Eq. (fT4"|) as 



Hid = I dz \ + huj r a s ip ip ijnp - ^hu r a 2 s ip t/j $ ^i/rip | . (16) 

1 2m Oz oz 2 / 



The first and second terms in Eqs. (|14l) or (|16p correspond to the Lieb-Liniger model. The third (cubic) term stems 
from the effective three-body collisions mediated by virtually excited states l . This third (cubic) term in Eq. \W\ 
violates the integrability in the ID system. 

The fact that the effective three-body interactions are dominated by virtual excitations of the lowest even-parity 
excited state may seem surprising, since the correct calculation of the effective two-body ID coupling constant requires 
taking into account of the infinite number of states [l6j]. However, in the latter case one deals with the calculation of 
the two-body wave function, which has in ID a 1/z singularity (stemming from the 3D boundary condition at r — > 
that provides the correct asymptotic form of the scattered s-wave). The removal of this singlularity yields the regular 
part of the two-body wave function and, through this regular part, the scattering amplitude and, hence, the effective 
coupling in ID. On the other hand, if one tries to obtain the effective ID interaction by adiabatic elimination of all 
the excited states, one gets a divergent series X)^lo n_1 ^ 2 m ^ ne expression for the effective ID coupling constant. 
The aforementioned regularization of the wave functions formally corresponds to the renormalization of this divergent 
series via substituting it by a finite expression lim^oo (J2n=o n ^ 2 — Jq dv i^ 1 / 2 ) [lj|. In our case, processes related 
to three-body collisions do not give rise to additional singularities in the many-body wave function, and no additional 
regularization is needed. The sum over all excited states thus converges. The convergence is rapid enough to ensure 
fair estimation of the whole sum by its first term. 

Accurate calculation of the effective three-body interaction potential Uzxs yields 

U3b(zi,z 2 ,z 3 ) = --fujJ r o^[y(z 1 ,z 2 ;Z3) +y(z 1 ,z 3 ;z 2 ) + y(z 2 , z 3 ; Zi)], (17) 
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where 



y( Zl ,z 2 ; z 3 ) = [6(2* - z t ) + 6(z 3 - z 2 )} £ ^ exp X A . (18) 

p—l * P r \ T J 

Obviously, the sum in Eq. (|18p converges and gives a sharp-peaked function rapidly (exponentially) decreasing at 
distances much larger than l r . Since Eq. ((T|) holds, we consider scattering events with transferred momenta much 
less than fi/l r . In this case we can use approximation y(z\, z 2 ; z 3 ) s» ^[d(z 3 — zi) + 6(z 3 — z 2 )]6(zi — z 2 ). Then, by 
taking the matrix element of the effective interaction Usb and dividing it by 3! (the number of permutation of three 
identical particles) we obtain the last term in the second-quantized Hamiltonian (|14|) . 

Before continuing we want to point out similarities with other recent works: (1) The mechanism discussed here 
is to a certain extent similar to the virtual association of atoms to a molecular dimer [l8j . In our discussion here, 
virtual excitation of radial modes during a two-atom collision temporarily localize the interatomic distance on the 
length scale ~ l r . Scattering a third atom on such a transient structure of finite size and mass 2m leads to an effective 
threebody collision. In Ref. [l8j], collisions of a third atom bring "virtual" dimers, enhanced in size by a Feshbach 
resonance down to the energy shell, thus bringing about "quantum chemistry" in ID. (2) In a similar way effective 
three-body interactions between polar molecules emerge, due to virtual transitions to an off-resonant internal state 
191. 



1 Note that in our previous work Il4[ the coefficient in front of this cubic term was estimated by a factor 4 too large. 
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B. Variational approach 



The cubic term in Eq. (| 16[) is negative and thus supports no bound ground state. Therefore we have to consider it 
as a first correction term to the purely pairwise interaction energy in the effective ID Hamiltonian. Considering the 

mean-field limit of Eq. (fT6|) . ip rj Jn\r> exp(i0), we get the energy density (per unit length) 



2m 



dz 
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Tz 



huj r a s nl D - ^fiiu r a 2 s n\ D , 



(19) 



this expansion is correct in the limit of the small linear density nma s < 1. If we neglect the contribution of the 
radial levels with the main quantum number larger than 2 by setting £ ss 1, we see that Eq. (|19[) is the expansion up 
to the cubic term of the energy density obtained by the variational method by Salasnich, Parola and Reatto [20| 
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Here one assumes the wave function of the transversal atomic motion to be oc exp[— (x 2 + y 2 )/(2a 2 )}, with a 2 = 
—^—y/1 + 2nnja s which minimizes £ var . 



IV. CALCULATIONS OF THE COLLISION RATES 



We now turn to the collision rates for these effective three body collisions. We start with calculating the rate Tk 1 k 2 k 3 
for the decay of a specific state \ki, k%, k 3 ) = a^o^ a£ |vac), |vac) being the vacuum state of the atomic field wherein 
atoms are absent (should not be confused with the vacuum of elementary excitations) due to three-body collisions. 
The final states of the decay are Ifc^fe^fejs) = a k' a k' a k' l vac )- To make the calcualtion simple we assume the ID 
bosonic gas being strongly non-degenerate (fc^T much higher than the chemical potential) and weakly-interacting. 
This enables us to neglect the probability double (and higher) occupation of any fc-mode, therefore assuming all 
involved atomic momenta to be different, and assume the elementary excitations coinciding with the atomic plane 
waves with the free-particle (quadratic) dispersion law. Then Fermi's golden rule yields 
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(21) 



To account for the condition k[ + k 2 + k' 3 — k\ + k 2 + k 3 we add an additional integration over k[ with the necessary 
delta-function: 



kxk 2 k 3 



L 2 
2irh 



{{{ dk[dk' 2 dk' 3 S 
J J Jw> 



2L 2 " 




(k[,k' 2 ,k 3 \ ^2 al 1 al 2 al 3 a qi a q2 a q3 \ki,k2,k 3 ) 



(22) 



Since the atoms are indistinguishable, we need to integrate over the volume W' in the fc'-space that corresponds 
to a unique ordering of the variables, e.g., k[ > k' 2 > k' 3 . The integral over W' of a function fully symmetric over 
permutations of ki,k 2 , k 3 amounts to 1/3! of the integral over the whole /c'-space. If all the involved momenta are 
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different (as is the case for a system far from degeneracy) all possible ways of ordering three bosonic creation operators 
with lower indices k[ , k' 2 , k' 3 and three bosonic annihilation operators with lower indices k\ , k 2 , k 3 finally yield 

{k[,k' 2 ,k 3 \ ^al 1 al 2 al 3 a qi a q2 a q:i \k 1 ,k2,k3) = (3!) 2 (23) 

To evaluate the integral in Eq. (|23p. we perform an orthogonal transformation from k^, k 2 , k 3 to Jacobi co-ordinates 
(here we deal with the ID analog of the hyperspherical co-ordinates, which are used in the three-body problem [21() 
in the wavenumber space: 



K — 7yg(^l + ^2 + ^3), ^12 — (^i ~ k 2 ), ^321 ~ \^ f^£S ~ 



(24) 



Further we introduce the hyperangle x' y i a 

k[ 2 = fc'sinx', k' 321 = fc'cosx', (25) 

and express Eq. (|2"3")l as 

6(V3k' c -ki-k 2 - k 3 )S(k' c 2 + k' 2 -k\- k\ - k\) = i\C 3b j^, (26) 

and accurate evaluation of the integral yields (cf. [lj|) 

C 3b = 3V3^ 2 ~ 6.88. (27) 

To calculate the three-body collision rate T 3 b per atom, we need to multiply Tk^H by the product of populations 
Nfkj of the states \kj) = a|..|vac), j = 1,2,3 (the occupation probabilities are normalized to unity, dk /& = 1) 
and integrate the whole k\,k 2 , fc3-space, divide by 3! to take into account the indistinguishability of the bosons. Then 
we obtain the number of three-body collisions per unit time in the whole system. This number should be divided by 
N to obtain the rate per atom: 



r 



36 



C 3 b -feV^ = C 3b C 2 ^. (28) 
m \ If J 



The result of Eq. (|28[) may seem counterintuitive at first: the collision rate is independent of temperature, and it is 
proportional to the dimensional parameter Q 2 and the radial confinement ui r . 

The physics behind the first observation is related to the fact that the collision kinetic energy is small compared 
to the virtual excitation energy. This was one of our assumptions in deriving the effective three body collisions and 
is required by the condition {ksT < hui r ) to be fulfilled when building a ID system (Eq. ([1])). Consequently the 
composite matrix element of the second-order process should not depend in leading order on the velocities (energies) 
of the colliding particles and hence on temperature (see Eq. ([9])). In addition the phase space volume for the scattered 
particles is independent on the incoming momenta k\ , k 2 , and k 3 . 

The other terms can be motivated the following basic physics considerations: Since effective three-body elastic 
scattering is the dominant process the scattering rate must be proportional to the 3D density squared: (nip/l 2 ) 2 . 
Furthermore, the scattering rate contains the square of the matrix element corresponding to the diagram in Fig. 1(a), 
where each vertex is proportional to a s , therefore this rate is proportional to af. The factor Ti/m provides the correct 
dimensionality (s _1 ). 

We can now compare the scattering rates for the two-body collisions T 2 b (Eq. (J7J) or effective three-body collisions 
^3b (Eq. (|28p) and evaluate their contributions to thermalization and the breakdown of integrability in ID systems. 
For ksT < TiLo r we find a simple scaling: 

r 3b C 3b 3\/3£ 2 

— = — (e k B T = (e kBT ~ 2.43 (e^ . (29) 

The relative importance of two-body collisions and the effective three-body collions mediated by virtual excitations is 
determined by the dimensionless quantity C,e k B T . For large C and small temperatures (fc^T <C hu> r ) the three-body 
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FIG. 2: Ratio between the scattering rates for the two routes to thermalization in quasi-lD systems: for the effctive three- 
body collisions and r 2 b for two-body collisions lea ding to excited transverse states. The points represent the ratios evaluated 
for various sets of experimental parameters from [ill ] (points), [l2( | (crosses), p]§| (triangle), and (diamonds). In these 
experiments the parameter £ was often close to 0.007 (the ratio T3b/T2b for £ = 0.007 exactly is shown by the solid curve). For 
comparison, we plot also T 3b /T 2 b for £ = 0.002 (dashed curve) and 0.02 (dot-dashed curve). Units on the axes are dimensionless. 



scattering rate due to virtual excitations dominates, and can lead to thermalization even when the thermalization 
processes due to two-body collisions are frozen out. For example, in typical atom chip experiment [lll - tl3T ] the three 
body reate T^i, dominates above the two-body collisions at fc^T < ~Kw r . A detailed comparison of the two rates F2& 
and and their relation to typical experimental parameters is given in Fig. [3] The scattering rate due to virtual 
excitations of the radial modes can dominate over real excitations for typical parameters of the recent experiment 

a 

The above calculation was for a non-degenerate ultracold gas. In a degenerate gas we need to consider the 
Bogoliubov-type spectrum of elementary excitations [22], that is phononic in the long- wavelength regime, as well 
as the relation between atoms and elementary excitations via the Bogoliubov transformations and bosonic amplifica- 
tion of scattering to modes, which are initially occupied. 

Taking into account all these factors, we find the rate of damping of a fast particle in a quasicondensate (see also 
Ref. [13): 

r da mp = 9 ^VC 2 J 00 dK £ dx , j°° d k>k'(l + n K )(l+n k ,)(l + n k ,)x 

S^Sk^Sk'Jik^ - ko/VfyS^ + T] k ' 2 + r)k> 3 - %„)• (30) 
Here the momenta of the scattered elementary excitations are defined by the expressions reciprocal to Eq. (IM1 l2~5t : 

k{ = -^ + y|fc'cos( X '-2V3), ^ = ^| + y|fc'cos(x' + 27r/3), 

k' 3 = -^ + y|fc'cos X '. (31) 
The energy of a mode with the momentum fik is Ek = h 2 rjk/(2m) with 



% = y/k*(k 2 +8n 1D a a /l*). 
The static structure factor of a quasicondensate is 



Sk = k 2 /rj k . 
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In equilibrium the population of the elementary mode with the momentum hk is given by the Bose-Einstein statistics 
with the mean occupation number for the mode with the momentum hk 



1 

nk 



exp[e fe /(fc B T)] -1' 

The initial kinetic energy % 'fc 2 , / '(2m) of the fast atom is assumed here to be large compared to both the mean- field 
interaction energy per particle in the quasicondensate and the temperature: 

k 2 > n lD a s /l 2 rl k 2 > mk B T/h 2 (32) 

Under condition (|32|) one of the scattered particles is always fast, and one of the three structure factors appearing in 
Eq. (|3U)l is always very close to 1 (and the corresponding occupation number is close to 0). In the most scattering 
events the other two particles are also fast, and, hence, Su ~ 1 and nu ~ for all three particles, j = 1,2,3. Only 
for the scattering events with small transferred momentum two of the structure factors are significantly less than 
1 and/or the corresponding populations approach the high-temperature limit ksT/ek'.- However, in the practically 
interesting case where k B T ~ huj r riiDa s the contribution of scattering events with small transferred momentum is 
relatively small, and 

rl: mp «^VC 2 . (33) 

The result of Eq. (|33|) for a fast atom in a quasicondensate is also obtained by Tan, Pustilnik and Glazman [24| . 



V. CALCULATIONS OF THE THERMALIZATION RATES 



We now turn to quantify thermalization in tightly confined Bosons in a quasi-lD geometry by both two body 
collisions and the effective interaction (|14[) . We again for simplicity consider a non-degenerate, weakly-interacting 
(the Lieb-Liniger parameter [5] 7 = 2a s / (riml 2 ) being much less than 1) gas of bosonic atoms. The assumptions of 
weak interaction and non-degeneracy enable us to express the three-particle distribution function through the product 
of single-particle distribution functions In contrast, calculation of relaxation via three-body collisions of low-energy 
excited states in the stronger interacting regime and especially for 7 > 1, requires to consider the (strong) quantum 
correlations in the quasi-lD bosonic system. We will discuss the effects of correlations on scattering rate and on 
thermalization in section IVII 

We start be writing the Boltzmann equation with a three-body collision integral [25j . taking into account the 
indistinguishability of the particles: 

Jt fk = T 36 J dk 1 J dk" J £ (f Ko fK_J K+1 - fkfvfw) , (34) 

with 

K s = k + k + k + y|fccos(7 + 2 S7 r/3), s = 0,±l, (35) 

k = ^ fc 2 + k' 2 + k" 2 - (fc + ^ ^ . (36) 

Eq. (|34l) can be easily understood: After integration over k' and k" the loss term in Eq. (|34"|) is simply — T^fk, which 
is the elastic three-body collision rate per atom. On the other hand, the three-atom state \k,k',k") is populated 
by elastic three-body collisions from those states \Kq, if— 1, K+x) which have the same center-of-mass momentum: 
Kq + if_j + if + x = k + k' + fc". Since the kinetic energy of the relative motion, % 2 k 2 / (2m) is conserved, the states 
|ifo, if— 1, if+i) (from where the state |fc, fc', fc") can be populated from) can be fully parametrized by the hyperanglc 
7- 

We now use the following Ansatz for the perturbed momentum distribution 

f k {t) = ^l^eM-k 2 /k 2 h )ll + e 4 (t)H 4 (k/k th )}, (37) 



10 





12 


80/ 70 / / 




10 


/ / 60 / / 


7 


8 




-O r- 1 
CO 21 

U 


6 






4 


^^^^^^^^^^^^^^^^^^ 




2 





3 4 5 6 7 

U) r l2n [kHz] 



FIG. 3: Dependence of the rate Tf%, of thermalization induced by effective three-body collisions in a weakly-interacting, quasi- 
1D 87 Rb gas on the radial trapping frequency for the linear densities hid from 80 fj,m to 30 fj,m (from top to bottom) with 
the step 10 /mi -1 (dot-dashed curve). 



with the numerical constant 



to solve Eq. Therby k^ — y / 2mkBT /% and H4 is the Hermite polynomial of the 4th order. This Ansatz is the 

simplest nontrivial perturbation that retains J dk kfk = 0. We then proceed to linearizing Eq. (|34p with respect to 
the perturbation amplitude £4,(t) and obtain an exponential solution e^{t) — £4(0) exp(— rj-Jt) with 

r W = C ^{^f) 2 < = C l^e (38) 

q 4] = lc„ = ^« 2.04. (39) 

To estimate the validity of our Ansatz we note that using a higher-order Hermite polynomial H n in Eq. (|37[) 
leaving the functional dependence on the parameters of the system unchanged and leads only to a minor modification 
of the numerical prefactor. For example, for n = 5 and 6 the thermalization rates are given by ^63;, w r C, 2 and 
lyC^u;,, C 2 , respectively. It is interesting to note that in these ID systems the thermalization rate due to the effective 
three-body collisions (TS) is about a factor 3 smaller then the collision rate (Tab). This suggests that in ID systems 
thermalization requires also about 3 collisions, similar to 3D (26[. Fig. |3] shows numerical values of r?£ as a function 

of the ID density of 87 Rb atoms and the radial trapping frequency. 

For comparison we calculate numerically the thermalization rate T 2 ^ for two-body collisions involving the real 
transitions between the ground and excited radial states. We follow hereby the same Ansatz and perturb the velocity 
distribution of atoms in the ground and excited state as given by Eq. (|37l) , the Boltzmannian distribution of overall 
populations between the levels being kept intact. In the parameter range of interest we find numerically rj^ 
(0.33 ±0.03)r 2& , i.e. 

2huj r 

« 0.93w r Ce~~. (40) 

The ratio of the thermalization rates for the two-body and three-body processes is therefore very close to the respective 
ratio of the collision rates, shown in Fig. [5] 

It is interesting to note that we find for both processes that thermalization in ID needs about 3 collisions capable 
to distribute energy. This is very close to the 2.7 collisions required for thermalization in 3D (26j . 

For the typical parameters of an ultracold 87 Rb gas on an atom chip [lj| (oj 7 . w 2tt x 3 kHz, nm ~ 50 /im -1 ) we 
obtain rj^ ~ 2 s _1 . This thermalization rate is temperature-independent and much larger than the one calculated 



11 



from the simple two-body collisions with the energy sufficient to excite radial modes r?£ «3x 10~ 3 s^ 1 at the lowest 
temperatures measured (30 nK). The estimated r 3 4 & j is consistent with the time needed for evaporative cooling of a 
87 Rb gas on an atom chip well below fwj r [HI, E2] • 



VI. SUPPRESSION OF THERM ALIZATION BY ATOMIC CORRELATIONS 

The thermalization rate r 3 4 & , given by Eq. (|38| was calculated for a weakly-interacting, non-degenerate gas. Cal- 
culation the thermalization rate r 3 4 j G in a general case requires to take into account additional physics. First we 
need to consider the effects of quantum degeneracy and second the fact that the disperison relations for the elemen- 
tary excitations in a ID quantum (degenerate) system may differ significantly compared to a free particle, especially 
for phonon-like excitations. These effects, together with the bosonic amplification of the scattering into thermally 
populated modes, tend to accelerate thermalization. A third observation is that the three body rates require three 
particles to be close to the same location. Quantum mechanically this is characterized by the third order correlation 
function gj, (0) . A full consideration of the above competing effects will require extended numerical analysis of many 
particular cases and transcends beyond the scope of the present manuscript. We will give here physical arguments of 
what to expect. 

We start by pointing out that the form of the secondary-quantized Hamiltonian Eq. (|16[) allows us to give a simple 
estimate of the ratio of these two rates 

rf 4 f/rf| = ^ 3 (o)/6) 2 . 

Here g is a phase-space factor accounting for the dispersion law of elementary excitation, which changes from free- 
particle-like to phonon-like. But it changes the thermalization rate less dramatically than the second factor associated 
with the local three-body correlation function 

.93(0) = (iP 1D (z)^ 1D {z))/nl D , 

and will become the dominating factor when approaching the strongly correlated regime (7 > 1). For a non-degenerate 
weakly-interacting Bose gas 93(0) = 3! = 6. For a degenerate ID Bose gas g% (0) has been recently calculated for the 
whole range of atomic repulsion strength (0 < 7 < 00) by Cheianov et al. [231. I n the zero-temperature limit 33(0) 
rapidly decreases from 1 to 167r 6 /(157 6 ) as 7 grows from to values 7 ^> 1. 

We now turn to the above conjecture on suppression of the thermalization by atomic correlations. A detailed 
calculation can be found in Ref. [15|, here we sketch the basic physics argument. To look at the correlations we start 
by considering N identical bosons in ID configuration with the Hamiltonian ITS] which, after rescaling of units, takes 
the form 

N d 2 

j=i 3 j>j> j>j'>j" 

c = 2a s /l 2 is the strength of interaction of two atoms in the tight waveguide with ground state size l r . Uyo is obtained 
by by adiabatic elimination of transverse modes virtually excited by the 3D short-range pairwise atomic interaction 
fl4j . The explicit form of [/3b is given by Eq. (p~7l) within a numerical prefactor, U-sx, = [Ti 2 / (2m)]Uzb- 

We follow now our detailed calcultions in Ref. [15[ and estimate the three-body scattering amplitude in the presence 
of the delta-functional pairwise repulsive interactions. The stronger the pairwise interparticle repulsion, the smaller 
is the probability of a close encounter of three particles which will result in a suppression of the three-body scattering 
amplitude. The simplest case is to analyze the Hamiltonian (14"TT) for N — 3 particles. For that purpose we express it 
in hyperspherical coordinates R, \ defined as [2l[ 

Z 1 +Z 2 +Z 3 z Y -z 2 [2 ( z 1 +z 2 \ . . 

Z c = — , Rsmx= ~ , Rco8x=y^\Z3 ^ — J (42) 



and obtain for the Hamiltonian 

3 

dZ 2 RdR^dR Wdx 2 ' If 



r) 2 1 r) r) 1 8 2 \F)r . . 

H - -7^-^ R 7^-W7^ + nrzZ^-^) + UMx). (43) 



u=-2 
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The corresponding Schrodinger equation for the three-particle wave function is 

H^(z!,z 2 , z 3 ) = (fcf + k\ + kl)^(z!,z 2 , z 3 ) ■ 

The wavenumbers kj are defined from the set of transcendental equations [5j , provided that the periodic boundary 
conditions are set on the interval of the length L. By setting L — > oo, we obtain a continuous spectrum, where kj's 
are real for repulsive interaction (c > 0). We can now separate the center-of-mass motion and describe the relative 
motion in hyperspherical coordinates. This leads to the Ansatz: 

$>(zi,z 2 ,z 3 ) = exp[i(fei + k 2 + k 3 )(zi + z 2 + z 3 )/3]ip T (R,x) 

where the kinetic energy of the relative motion is given by 

k 2 = i[(fci - k 2 f + (k 2 - fc 3 ) 2 + (k 3 - fci) 2 ] 

. In the adiabatic hyperspherical approximation [30j , which holds in the long- wavelength limit 

k < c (44) 

and neglects coupling of different scattering channels as well as accumulation of phase shifts of the scattered wave 
due to non-adiabatic effects we get 

tP t (R,x)=F (R)B (R,x)- (45) 

The hyperangular part B (R, x) of this wave function is the eigenfunction of the Hamiltonian (|43[) with fixed R, 
corresponding to the lowest eigenvalue \$(R), which is the smallest positive root of the transcendental equation 

X(R) tan[7rA(i?) /6] = cR/V2. (46) 

Using the regular hexagon symmetry group [3l[ of the Hamiltonian (|4"3")h we write the hyperangular part of Eq. (|4"5)l 

as 

BoCR, x) = B {R, x) + B {R, X - 2tt/3) + B (R, X + 2tt/3), (47) 

where 

B»Wri^{r |Ao(S)( ' /6 " lxl,, 'Ll™t W 

After integrating out the hyperangular variable, the Schrodinger equation in the adiabatic hyperspherical approxima- 
tion reduces to 



RdR dR 



rA ° (i?) ' Uoo(R) 



R 2 



F = k 2 F (49) 



where 



_ r 3 ^fo'wx)t/»wx) 

with the boundary conditions requiring Fq to be finite for both R = and R — > oo, for the "partial wave" corresponding 
to the lowest eigenvalue Ao(-R), whose asymptotic expressions are 



A ° (i?) "i v "S cR<<1 ■ ^ 

We can solving now Eq. (|4"9")l analytically in two regions, cR -C 1 and cR ^> 1, with Ao(i?) approximated by Eq. 
(f5Tj) . and tailoring the solutions by quasiclassical expressions for Fq(R) in the intermediate hyperradius range. The 
scattering amplitude /o (for its definition in planar geometry see [32l - [35j ) can then be obtained from the asymptotic 
form of the wave function at R — > oo 

F (R) w J 3 (kR) - ifoH^ikR), (52) 
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where H^\z) — J^{z) + iY^(z) is the Hankel function of the first kind and Jz(z), Y^(z) are the third-order Bessel 
functions. The behavior of F at R — > is defined by the details of the potential Uoo{R), but the result can be finally 
expressed via the effective vertex of the three-body elastic collisions, thus yielding 

A . (53) 

where 17 w 1 is a numerical constant and £ ss 1.15 is defined by Eq. (fTS"j) . We use thereby the fact that Uoq(R) differs 
significantly from zero on the length scale l r , over which a virtually excited particle can propagate, and cl r <C 1. 

From Eq. (|53|) we conclude that the three-body scattering amplitude decreases in proportion to (k/c) 6 as k/c — > 0, 
i.e., when the interaction is strong enough to induce significant atomic correlations. The three-body scattering rate 
in a ID system of bosons in the case of strong pairwise interaction is suppressed by a factor ~ (k/c) 12 . By averaging 
over collision momenta in a moderately-excited strongly-interacting state we obtain the scattering rate suppression 
factor- ((fc/c) 12 ) - 7 - 12 . 

We can now compare our result with the zero-distance three-particle correlation function 33(0). In the strong 
interaction limit 7 » 1 53(0) cx 7~ 6 [U and gives a direct physical motivation of our original conjecture Q that 
the pairwise interactions and the quantum correlations induced by them in a strongly-interacting ID bosonic system 
suppress the three-body elastic scattering rate, and, hence, thermalization, by a factor cx g!(0). In other words, strong 
quantum correlations extend the time scale, on which a quasi-lD system approaching the Tonks-Girardeau regime 
can be considered as approximately integrable. 

Moreover, we can corroborate this conjecture by observing that the thermalization rate is proportional to the square 

of the matrix element of the transition operator that is proportional to J dz ip ip ip ipipip. Since the energies of the 
products of the elastic three-body process are low (of about fceT), we may assume that the correlations in the initial 
and the final states are the same, and the transition matrix element can be regarded as proportional to (73(0) that 
yields again the (0) scaling of the rate. In contrast to this, the inelastic three-body processes are accompanied 
by a large energy release, and after an inelastic collision the newly formed dimer molecule and the fast atom leave 
the system almost immediately. Therefore the inelastic three-body relaxation rate in a ID ultracold Bose gas is 
proportional to the first power of 53(0) (37j . 

We now compare the calculated thermalization rates to the quantum Newton's cradle experiment [Toj where a 
lower boundary for the damping time of the ID motion towards a Gaussian profile was estimated. For three different 
Lieb-Liniger parameters 7 = 1.4, 3.2, 18 Kinoshita et al estimate lower bounds to the thermalization time from the 
consistency of the the experimental momentum distributions with the experimentally observed heating during the 
time interval of 0.5, 0.5, 1.0 seconds probed. They find one sigma lower limits of 2.6 s, 25 s, and 13 s respectively. 
In their experiment the motion of two groups of 87 Rb atoms was excited at the relative velocity equal to 4 recoil 
velocities, which is far above the width of the ground-state velocity distribution. We therefore can not expect the 
collision rate to be suppressed in proportion to g 2 (0) that, as described above for slow collisions. Instead, we have to 
apply the estimates for damping of a fast particle discussed in the end of section IV. In a strongly-correlated system 
one has to take the particle correlations into account and eq. [33] has to be multiplied by a factor 32(0) denoting the 
two-particle correlation function at zero distance (Tan, Pustilnik and Glazman [241]): 

C^T^C'ftfO). (54) 

Substituting the experimental parameters of Ref. [13] and taking the values for 52(0) from Ref. [38|], we obtain 
r fc( ^ mp s» 15 s _1 , 1.7 s _1 , and 7 x 10~ 3 s _1 for 7 = 1.4, 3.2, and 18, respectively. To compare this calculated damping- 
rates to the experiment in [Io| one has to consider that (1) the two colliding clouds overlap only for a very short 
time during each oscillation and (2) that the thermalization rate is a factor 3 longer (section V). Taking this into 
account we estimate the respective thermalization times of 2.6 s, 35 s and > 1000 s. These rates are consistent with 
the experimental findings of Ref. [Io| . For a more detailed comparison one would need longer time scale experiments 
with lower intrinsic heating and more detailed calculations of the dynamics including the damping due to three body 
collisions discussed here and in 12411 . 



VII. CONCLUSION 



A radially confined atomic gas is never perfectly ID, and radial motion can be excited, either in reality or virtually 
even if Eq. ([T]) holds. This possibility leads to effective three-body collisions, which arise in the second order of 
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perturbation theory and can be associated with virtual excitation of radial modes. These processes lead to thermal- 
ization even when two body collisions are frozen out at ksT <C Tiuj r and provide a mechanism to break integrability 
in ID systems. In other words, the freeze-out of the radial modes is only a necessary, but not sufficient condition for 
integrability in ID systems. Our estimations of the relaxation rates for weakly interacting; quasi-lD Bose gases are 
consistent with recent experimental observations for weakly-interacting quasicondensates [ill [l2| . 



These effective three-body collisions can be suppressed by quantum correlations caused by strong pairwisc repulsions. 
If they dominate, as in a strongly correlated ID Tonks- Girardeau gas, they suppress the influence of the integrability- 
breaking interaction term. The thermalization rate decreases in proportion to g§(0) as the system enters the regime 
of strong correlations (7 > 1), and the system (remaining non-integrable in the strict sense) behaves like (almost) 
integrable on time scales short compared to the inverse thermalization rate. 

The effective three-body collisions, and their suppression by quantum correlations should be accessible in experi- 
ments looking at the damping of fast, particle- like excitations in systems with 7 varying in a broad range of values 
from less than 1 to ~ 10. 
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